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ABSTRACT 

The  unique  nature  of  the  molecular  structures  of  gas  hydrates  results  in  curious  acous¬ 
tic  properties  which  have  yet  to  be  adequately  characterized.  Understanding  the  acous¬ 
tic  behavior  of  hydrates  in  liquids,  in  bubbly  liquids,  and  in  sediments  containing  liquids 
and/or  gas  is  vital  for  surveying  their  location  using  seismic  or  echosounding  techniques  and 
may  become  a  key  tool  for  monitoring  hydrate  dissociation  and  its  possible  link  to  climate 
change.  Acoustic  properties  of  gassy  substances  are  known  to  have  a  strong  dependence 
on  excitation  frequency;  however,  tabulated  values  of  hydrate  sound  speeds  are  most  often 
measured  at  high  frequencies  (>200  kHz)  despite  modern  location  methods  which  use  fre¬ 
quencies  below  100  kHz.  This  presentation  details  a  laboratory  experiment  in  which  the 
dissociation  pressures  of  natural  structure  I  and  structure  II  methane  hydrate  samples  were 
determined  by  measuring  their  low-frequency  acoustic  velocity  in  a  liquid  as  a  function  of 
hydrostatic  pressure.  [Work  supported  by  The  United  States  Office  of  Naval  Research.] 
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NOMENCLATURE 

c  speed  of  sound 

/  frequency 

G  shear  modulus 

K  bulk  modulus 

L  resonator  length 

n  mode  number 

P  hydrostatic  pressure 

S  salinity 

si  structure  I  gas  hydrate 

sll  structure  II  gas  hydrate 

T  temperature 

VS  A  vector  signal  analyzer 

A  acoustic  wavelength 

p  density 

X  volume  fraction 


INTRODUCTION 

Gas  hydrates  are  often  formed  and  found  in 
ocean  sediments  along  continental  margins.  The 
compounds  are  also  known  to  form  as  a  skin  on 
rising  methane  bubbles,  as  reported  by  Rehder 
et  al.  [1],  Heeschen  et  al.  [2],  and  Sauter  et  al.  [3]. 
To  improve  seismic  detection  and  characteriza¬ 
tion  of  hydrates,  it  is  necessary  to  gain  a  better 
understanding  of  their  acoustic  behavior  in  the 
variety  of  media  in  which  they  are  found.  Of  par¬ 
ticular  import  is  low-frequency  acoustic  behavior 
of  hydrates,  which  has  yet  to  see  extensive  inves¬ 
tigation.  As  a  step  toward  a  better  understand¬ 
ing  of  the  behavior  of  hydrates  in  acoustically 
complex  media  such  as  in  water-saturated  sed- 
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iments,  this  work  focuses  on  the  low-frequency 
acoustic  behavior  of  structure  I  and  structure  II 
methane  hydrates  in  liquids  and  in  bubbly  liq¬ 
uids. 

Definitions 

Gas  hydrates,  often  referred  to  as  clathrate  hy¬ 
drates ,  gas  clathrates ,  or  simply,  hydrates,  are 
crystalline  molecular  host  structures  which  en¬ 
trap  guest  molecules  without  chemical  bonds. 
The  substance  forms  in  high-pressure,  low- 
temperature  environments,  and  requires  such 
conditions  to  remain  stable.  The  multiple  names 
for  the  compound  are  often  used  interchange¬ 
ably;  however,  each  term  describes  specific  traits 
of  the  structure.  Powell  defined  clathrate,  from 
the  Latin,  clathratus,  meaning,  “to  encage,”  as 
a  general  term  for  any  such  cage-like  molecu¬ 
lar  structure  [4],  A  clathrate  hydrate  is  then  a 
compound  in  which  water  forms  the  rigid  lattice, 
and  gas  clathrate  hydrate  specifies  that  the  guest 
molecule  in  the  structure  is  a  gas.  In  this  pa¬ 
per,  the  term  clathrate  will  be  dropped  and  the 
compounds  will  be  referred  to  as  hydrates,  gas 
hydrates,  or  methane  hydrates,  dependent  upon 
the  appropriate  level  of  generality  or  specificity. 
Three  molecular  structures  of  gas  hydrates  are 
known  to  exist.  Structure  I  is  a  body-centered- 
cubic  orientation  which  forms  naturally  in  deep 
oceans  from  biogenic  gases.  It  is  worth  not¬ 
ing  that  this  molecular  geometry  can  trap  great 
quantities  of  gas  in  a  small  volume,  concentrat¬ 
ing  methane  by  a  factor  of  164  relative  to  its 
volume  at  STP  [5].  Structure  II  is  a  diamond 
lattice  within  a  cubic  framework,  formed  from 
thermogenic  gases,  often  in  oil-rich  areas  [6]. 
Structure  H  is  a  rare  hexagonal  form,  which  has 
only  recently  been  shown  to  exist  in  nature  [7]. 
The  focus  of  the  present  work  is  on  structure  I 
and  structure  II  hydrates  hosting  molecules  of 
methane  gas. 

DESCRIPTION  OF  EXPERIMENT 

A  pressure-  and  temperature-controlled  cham¬ 
ber  was  designed  to  investigate  the  effect  of 
hydrostatic  pressure  on  the  acoustic  properties 
of  brines  containing  samples  of  methane  hy¬ 
drates.  An  annular  cylindrical  borosilicate  res¬ 
onator  with  inner  diameter  of  52  mm,  outer 
diameter  of  70  mm,  and  length  459  mm,  was 


Figure  1:  Schematic  of  the  experimental  appa¬ 
ratus. 

wrapped  in  coiled  copper  tubing  through  which 
chilled  ethylene  glycol  at  a  temperature  of  — 20°C 
circulated.  The  exposed  copper  tubing  was  then 
covered  with  fiberglass  insulation.  The  bottom 
of  the  resonator  was  sealed  with  a  nitrile  rub¬ 
ber  membrane  and  the  resonator  rested  atop  an 
open-cell  foam  block,  creating  an  approximately 
pressure-release  boundary  condition  at  the  lower 
terminus.  An  audio  amplifier  powered  an  elec¬ 
tromagnetic  shaker,  to  which  an  aluminum  pis¬ 
ton  was  attached.  The  piston  was  placed  near 
the  top  of  the  column  of  brine  and  oriented  to  ex¬ 
cite  longitudinal  acoustic  modes  within  the  fluid 
enclosed  by  the  resonator  walls  while  a  small 
laboratory  hydrophone,  with  its  cable  encased 
in  a  water-filled  stainless  steel  sheath,  sensed 
the  acoustic  response  of  the  system.  The  hy¬ 
drophone  signal  was  amplified  and  bandpass  fil¬ 
tered  with  a  voltage  preamp.  A  vector  signal 
analyzer  (VSA)  generated  the  acoustic  signals 
(band-limited  periodic  chirps)  and  digitized  the 
acoustic  spectrum  of  the  sample  for  each  mea¬ 
surement.  The  apparatus  and  procedures  used 
in  the  present  work  closely  mimics  the  exper¬ 
imental  design  described  by  Wilson  et  al.  in 
Ref.  [8].  A  schematic  of  the  apparatus  is  given 
in  Fig.  1. 

A  copper  wire  busbar  cage,  shown  in  Fig.  2, 
was  constructed  to  suspend  six  approximately 
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Figure  2:  Schematic  of  borosilicate  resonator 
with  copper  wire  cage  holding  six  hydrate  sam¬ 
ples.  Open-cell  foam  is  seen  at  the  lower  termi¬ 
nus  of  the  resonator. 


equally-spaced  pieces  of  gas  hydrate  in  the  res¬ 
onator.  The  samples  float,  so  to  keep  them 
in  place,  the  cage  was  wrapped  in  copper  wire 
mesh.  To  prevent  bubbles  from  attaching  to 
the  copper  cage,  the  cage  was  coated  with  a 
hydroxy-terminated  polydimethylsiloxane  solu¬ 
tion  and  allowed  to  dry  before  it  was  placed  in 
the  resonator. 

Natural  samples  of  methane  hydrates  were  col¬ 
lected  through  coring  operations  between  1998 
and  2000  and  stored  in  liquid  nitrogen  at  at¬ 
mospheric  pressure  until  January  2010.  At  that 
time,  the  hydrates  were  packed  in  a  dewar  with 
liquid  nitrogen  and  shipped  from  the  storage  fa¬ 
cility  at  the  Naval  Research  Laboratory  in  Wash¬ 
ington  D.C.  to  Applied  Research  Laboratories  in 
Austin,  Texas.  Geographical  origins  and  molec¬ 
ular  structures  of  the  samples  are  listed  in  Ta¬ 
ble  1. 

For  each  test,  a  hydrate  sample  was  divided  into 
six  pieces,  each  approximately  2  cm  in  diameter, 
and  placed  in  the  copper  cage.  The  cage  was 
lowered  into  the  pre-chilled  resonator  which  was 
filled  with  a  brine  of  salinity  S  =  131  ±  2%o. 
With  the  resonator  in  place  inside  the  pressure 
vessel,  the  piston  of  the  shaker  and  hydrophone 
were  positioned  near  the  top  of  the  brine,  and 
the  vessel  was  quickly  sealed  and  pressurized 
with  nitrogen  to  approximately  2.5  MPa.  The 
height  of  the  liquid  column  inside  the  resonator, 
L,  was  determined  by  viewing  a  video  image 
from  a  camera  placed  inside  the  pressure  ves¬ 


sel  and  aimed  at  the  top  of  the  resonator.  The 
VS  A  recorded  an  acoustic  spectrum  given  by 
the  linear  frequency-domain  average  of  30  trans¬ 
fer  functions  (between  the  hydrophone  and  ex¬ 
citation  signals)  over  a  frequency  range  which 
was  determined  based  upon  the  resonance  fre¬ 
quencies  observed  for  each  sample  and  hydro¬ 
static  pressure.  After  recording  each  spectrum, 
a  small  amount  of  gas  was  vented  from  the  cham¬ 
ber  to  slightly  reduce  hydrostatic  pressure  in  the 
vessel.  Measurements  were  repeated  at  decre¬ 
menting  pressures  until  the  acoustic  signal  was 
no  longer  discernible.  Shifts  in  resonance  fre¬ 
quencies  were  often  observed  during  pressure  re¬ 
duction;  however,  no  spectral  changes  were  seen 
more  than  10  seconds  after  pressure  reduction. 
To  ensure  stability  of  the  system,  measurements 
were  performed  120  seconds  after  each  pressure 
reduction. 

For  each  recorded  spectrum,  sound  speed  of  the 
mixture  was  calculated  by  the  slope  method  de¬ 
scribed  in  the  following  section,  taking  only  the 
two  lowest-frequency  resonance  peaks.  At  high 
hydrostatic  pressures,  several  resonance  peaks 
were  distinct  over  the  frequency  range  and  a  lin¬ 
ear  fit  of  the  resonance  frequencies  versus  mode 
number  showed  a  low  coefficient  of  determina¬ 
tion.  However,  at  low  hydrostatic  pressures, 
higher-mode  resonances  faded  and  became  in¬ 
distinguishable.  Thus,  for  consistency,  only  the 
two  lowest  resonance  frequencies  were  used  for 
all  spectra. 

The  experimental  design  was  tested  by  measur¬ 
ing  the  sound  speed  of  a  brine  (S'  =  98%o) 
chilled  to  —  0.3°C.  Without  the  copper  wire  sam¬ 
ple  holder  in  the  resonator,  in-resonator  sound 
speed  was  measured  as  1347  m/s.  After  account¬ 
ing  for  the  elastic  waveguide  effect  of  the  borosil¬ 
icate  tube  (described  below) ,  the  freefield  acous¬ 
tic  velocity  of  the  brine  was  determined  to  be 
1546  m/s,  approximately  0.4%  greater  than  the 
predicted  value.  This  error  is  likely  caused  by 
uncertainties  in  the  measurement  of  L  or  tem¬ 
perature  in  the  column.  Next,  the  empty  copper 
cage  was  placed  in  the  resonator  with  the  brine 
and  the  measurement  was  repeated.  The  pres¬ 
ence  of  the  cage  in  the  resonator  had  minimal 
effects  on  the  overall  structure  of  the  acoustic 
spectrum  of  the  water  column,  but  resulted  in 
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sample 

structure 

origin 

T  [°C] 

CH4  [%] 

c2h6  [%] 

c3h8  [%] 

c4h10  [%] 

CMsI 

I 

Cascadia  Margin 

-8.3 

99.79 

0.21 

- 

- 

CMsII 

II 

Cascadia  Margin 

-4.7 

81.89 

10.4 

3.7 

1.40 

HMsI 

I 

Hakon  Mosby 

-6.2 

99.49 

0.14 

0.12 

0.15 

GMsII 

II 

Gulf  of  Mexico 

-6.5 

29.67 

15.35 

36.61 

13.61 

Table  1:  Structures,  geographical  origins,  temperatures  and  compositions  of  hydrate  samples.  Com¬ 
positions  were  determined  using  a  Varian  Saturn  2200  mass  spectrometer,  then  quantitatively  deter¬ 
mined  with  a  Varian  CP  3800  gas  chromatograph,  as  described  in  Ref.  [9].  Percentages  of  butane  and 
isobutane  have  been  combined  in  the  C4H10  column. 


a  0.3%  decrease  in  the  measured  sound  speed  of 
the  brine. 

EXPERIMENTAL  MODELS  AND 
DATA  ANALYSIS 

Measurements  performed  in  the  present  ex¬ 
periment  employed  an  apparatus  known  in 
the  acoustics  community  as  a  one-dimensional 
acoustic  resonator  or  waveguide.  Sound  speeds 
of  hydrate/brine  mixtures  measured  in  the 
resonator  were  numerically  corrected  to  ac¬ 
count  for  a  well-known  systematic  error  due 
do  elasticity  of  the  resonator  walls,  then  the 
data  was  further  analyzed  to  extract  material 
properties  of  the  hydrate  samples  using  Wood’s 
model  of  sound  propagation  through  multiphase 
mixtures.  Changes  in  sound  speeds  of  the 
mixtures  observed  as  a  function  of  changes  in 
hydrostatic  pressure  offer  clear  measurements  of 
phase  boundaries,  and  are  in  rough  agreement 
with  established  hydrate  stability  models. 

1-D  Resonator  Sound  Speed  Measurement 

The  low-frequency  sound  speed  of  a  fluid  can 
be  measured  by  determining  the  resonance  fre¬ 
quencies  of  the  fluid  sample  inside  an  acous¬ 
tic  waveguide.  This  technique  facilitates  lab¬ 
oratory  measurements  at  low  frequencies  be¬ 
cause  the  resonator  length  need  only  be  one- 
half  acoustic  wavelength,  as  compared  to  time- 
of-flight  measurements,  which  require  sample 
lengths  equivalent  to  several  wavelengths.  In  the 
present  model,  the  inner  walls  of  the  borosili- 
cate  tube  are  assumed  to  be  rigid  and  the  ends 
are  approximated  as  pressure-release  boundary 
conditions.  Temporarily  neglecting  the  elastic¬ 
ity  of  the  borosilicate  tube,  the  expression  for 
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the  eigenfrequencies  of  the  resonator,  fn  = 
where  n  is  the  mode  number  and  L  is  the  length 
of  the  resonator,  yields  an  expression  for  the 
sound  speed  of  the  fluid  mixture, 


where  Afn/ An  is  the  slope  of  the  line  given  by  a 
linear  least  squares  fit  of  the  measured  resonance 
frequencies  versus  mode  number.  Returning  now 
to  the  finite  elasticity  of  the  wall  material,  the 
effect  of  the  wall  elasticity  is  seen  as  a  reduction 
in  sound  speed  relative  to  that  which  would  be 
observed  in  an  unconfined  volume  of  the  same 
material,  often  called  the  “free  held.”  An  exact 
analytical  model  of  sound  propagation  through 
a  finite-thickness  elastic-walled,  fluid-filled  cylin¬ 
drical  tube  is  described  in  Ref.  [10].  The  present 
work  involves  two  types  of  measurements.  The 
first  is  acoustic  measurement  of  hydrate  dissocia¬ 
tion  pressures,  which  relies  on  relative  changes  in 
sound  speed  due  to  the  presence  of  gas  bubbles. 
For  determination  of  dissociation  pressures,  the 
apparent  sound  speeds  observed  inside  the  res¬ 
onator  were  adequate,  and  the  correction  for  the 
elastic  waveguide  effect  was  not  applied.  Hy¬ 
drate  bulk  moduli  were  also  measured  using  the 
low-frequency  acoustic  technique.  For  these  ma¬ 
terial  property  measurements,  it  was  necessary 
to  correct  for  the  elasticity  of  the  tube  walls  us¬ 
ing  the  procedure  described  in  Ref.  [11]. 

Wood’s  Model  of  Sound  Propagation 
Through  Multiphase  Mixtures 

Measurement  capabilities  of  an  acoustic  res¬ 
onator  extend  beyond  the  simple  determination 
of  the  speed  of  sound  in  a  single  fluid  medium — 
Material  properties  of  individual  constituents 
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in  a  multiphase  mixture  may  also  be  obtained 
through  the  use  of  Wood’s  model,  which  is  de¬ 
scribed  here.  A  fluid  medium  may  be  approxi¬ 
mated  as  an  acoustically-homogeneous  mixture 
if  all  particles  constituting  the  mixture  are  much 
smaller  in  dimension  than  an  acoustic  wave  trav¬ 
eling  through  the  mixture,  the  particles  are 
evenly  distributed  throughout  the  mixture,  and 
the  particles  number  3>  1  per  wavelength.  By 
this  definition,  a  two-phase  medium  consisting 
of  a  liquid  hosting  a  distribution  of  gas  bubbles 
or  solid  particles  may  be  approximated  as  a  bulk 
medium  with  an  effective  bulk  modulus  Km  and 
density  pm.  The  effective  properties  of  the  mix¬ 
ture  are  obtained  as  a  linear  combination  its  con¬ 
stituent  properties  by  their  volumes, 


of  the  guest  particles  in  the  mixture, 


Xi  = 

V! 

(7) 

Vl  +  V2  +  V3 

X2  = 

v2 

(8) 

Vi  +  v2  +  v3 

X3  = 

v3 

(9) 

V  +  v2  +  v3 ' 

The  bulk  modulus  of  the  mixture  is  then 

I<  =  KlK2K*  (W) 

m  XiK2Kz  +  X2K1Kz  +  xzK1K2'  1  > 

where  K\ ,  AT,  and  it  3  are  the  bulk  moduli  of  the 
liquid  host  and  the  gas  and  solid  guests,  respec¬ 
tively.  Likewise,  the  density  of  the  three-phase 
mixture  is 


Km 


AT  if  2 

X.K\  +  (1  -  X)K2 


and 

Pm  =  (1  ~x)pi  +XP2, 


(2) 

(3) 


where  K\  and  AT  are  the  bulk  moduli  of  the 
liquid  and  guest  particles,  respectively,  and  p± 
and  p2  are  the  densities  of  the  liquid  and  guest 
particles.  The  volume  fraction  x  °f  the  guest  is 
given  by 


V> 

Vi  +  V2  ’ 


(4) 


where  V2  and  V\  are  the  total  respective  vol¬ 
umes  of  the  guest  particles  and  host  liquid  in 
the  resonator.  Wood’s  model  [12]  of  sound  speed 
through  the  mixture  is  then  found  as  a  function 
of  the  bulk  modulus  and  density  of  the  effective 
medium, 


Cm  — 


Cm  — 


'Am 

Pm 


or, 


IK1K2[(1-x)Pi  +  XP2]  1 
XK\  +  (1  -  x)AT 


(5) 

(6) 


Similarly,  Wood’s  model  may  also  be  employed 
to  predict  the  sound  speed  of  three-phase  media 
such  as  the  mixtures  of  seawater,  gas,  and  solid 
hydrates  which  have  been  observed  in  nature. 
For  the  three-phase  case,  volume  fractions  \q, 
X2 ,  and  X3  are  used  describe  the  relative  volumes 


Pm  =  XI  Pi  +X2P2  +X3P3-  (11) 

The  expressions  for  volume  fractions,  mixture 
bulk  modulus,  and  mixture  density  then  lead 
to  an  expression  for  the  mean  velocity  of  an 
acoustic  wave  traveling  through  the  three-phase 
medium, 


Cm  — 


I  K^Ksjxipi  +  X2P2  +  Xsps]  1 

X 1  AT  AT  +  X2ATA3  +  X3ATAT 


(12) 


Given  this  model  and  knowledge  of  the  ma¬ 
terial  properties  of  any  two  of  the  three  con¬ 
stituents  of  the  mixture,  it  is  possible  to  then 
determine  an  unknown  material  property  of  the 
third  constituent.  For  example,  the  bulk  modu¬ 
lus  and  density  of  a  water  or  brine  host  is  well- 
established.  Assuming  a  dissociating  methane 
hydrate  releases  bubbles  of  methane  gas  into  the 
mixture,  tabulated  properties  of  methane  gas 
should  suffice  as  values  for  AT  and  p2.  A  sim¬ 
ple  laboratory  measurement  of  the  hydrate  sam¬ 
ple  density  then  leaves  one  unknown  value,  AT 3, 
which  may  be  determined  by  rearranging  Eq.  12 
using  measured  values  of  Cm. 


Hydrate  Stability  Models 

In  recent  years,  several  groups  have  attempted 
to  model  the  stability  regimes  of  gas  hydrates. 
Parrish  and  Prausnitz  (Ref.  [13])  used  experi¬ 
mental  data  to  fit  constants  to  a  model  based  on 
the  work  of  Van  der  Waals  and  Platteeuw  [14]. 
Parrish  and  Prausnitz  stated  that  their  model 
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model 

relation 

Parrish  and  Prausnitz,  1972 
Dickens  and  Quinby-Hunt,  1994 
Sloan,  1998 

Peltzer  and  Brewer,  2000 

ln(P/101.325)  =  -1212.2  +  44344.0/T  +  187.719  In  T 
±  =  3.79  x  10~3  -  2.83  x  10”4[log10(P/1000)] 

P  =  exp(38.98  -  8533.8 /T) 

}  =  3.83  x  10-3  -  4.09  x  10-4[log10(P/1000)]... 
+8.64  x  10_5[log10(P/1000)]2 

Table  2:  Hydrate  stability  models  taken  from  the  literature.  Pressures  are  given  in  kilopascals  and 
temperature  is  in  units  kelvin. 


should  only  apply  to  structure  I  methane  hy¬ 
drates,  and  the  model  was  validated  for  a  tem¬ 
perature  range  of  0°C  to  27°C.  Two  decades  later, 
Dickens  and  Quinby-Hunt  developed  an  empiri¬ 
cal  equation  based  on  data  from  their  own  exper¬ 
iments  on  methane  hydrates  in  seawater  of  salin¬ 
ity  S  =  33.5%o  [15].  Sloan  later  gave  a  detailed 
explanation  of  the  thermodynamic  and  kinetic 
properties  of  each  of  the  three  known  structures 
of  methane  hydrates.  In  that  paper  he  gives  an 
expression  for  the  temperature-dependent  dis¬ 
sociation  pressure  of  methane  hydrate,  but  ne¬ 
glects  to  discuss  how  the  relation  was  developed 
or  which  hydrate  structure  it  describes  [6].  Soon 
after,  Peltzer  and  Brewer  fit  the  data  presented 
in  Sloan’s  paper  to  a  simple  second-order  polyno¬ 
mial  to  build  an  expression  which  was  modeled 
after  Dickens  and  Quinby-Hunt’s  relation  [16]. 
The  models  are  provided  for  reference  in  Table  2. 

The  four  models  discussed  in  this  paper  take 
various  forms,  but  each  describes  an  inverse  re¬ 
lationship  between  the  pressure  and  tempera¬ 
ture  of  hydrate  dissociation.  Of  these  mod¬ 
els,  only  the  Parrish  and  Prausnitz  model  is 
structure-specific.  Although  some  of  the  authors 
who  developed  the  models  listed  in  Table  2  ac¬ 
knowledged  that  the  dissociation  temperature- 
pressure  relation  will  be  affected  when  the  hy¬ 
drate  is  surrounded  by  seawater,  only  Dickens 
and  Quinby-Hunt  explicitly  stated  the  salinity 
of  the  host  liquid  considered  in  the  development 
of  their  models. 

RESULTS 

During  the  experiments  on  gas  hydrate  samples 
it  was  observed  that  the  acoustic  spectra  cap¬ 
tured  at  hydrostatic  pressures  above  500  kPa  to 
1000  kPa  showed  impressive  signal-to-noise  ra- 
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Figure  3:  Measured  low-frequency  sound  speeds 
of  hydrate  mixtures. 


tios  with  distinct  peaks  at  resonance  frequen¬ 
cies  and  deep  valleys  at  antiresonances.  Regu¬ 
larity  of  the  intervals  between  resonances  indi¬ 
cated  a  frequency-independent  sound  speed  be¬ 
low  9  kHz,  and  a  >50  dB  difference  between 
peak  and  valley  receiver  voltages  implied  a  low- 
noise  system  with  little  acoustic  loss.  As  hy¬ 
drostatic  pressure  in  the  chamber  decreased,  the 
hydrates  began  to  dissociate,  gas  formed  on  the 
samples,  and  the  bubbles  were  released  into  the 
brine.  The  introduction  of  bubbles  then  lowered 
the  effective  sound  speed  of  the  bulk  medium  in 
the  resonator,  thus,  shifting  resonances  lower  in 
frequency.  The  rising  bubbles  in  the  column  re¬ 
duced  signal-to-noise  ratios,  due  in  part  to  flow- 
generated  noise,  and  due  in  part  to  the  attenu¬ 
ation  of  acoustic  energy  by  conversion  into  heat 
during  oscillation  of  the  excited  bubbles.  Mea¬ 
surements  were  discontinued  when  the  acoustic 
signal  was  no  longer  discernible.  The  full  data 
set  and  a  detailed  description  of  the  experiment 
are  presented  in  Ref.  [17]. 

Measured  values  of  in-resonator  acoustic  veloc¬ 
ities  for  four  brine  mixtures  containing  the  gas 
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Figure  4:  Acoustically-measured  gas  hydrate  dissociation  pressures  compared  to  established  models. 
The  Parrish  and  Prausnitz  model  is  only  valid  for  comparison  to  results  obtained  for  structure  I  hy¬ 
drates.  Most  predicted  values  trend  lower  than  dissociation  pressures  measured  acoustically,  possibly 
due  to  the  models’  neglect  of  salinity. 


hydrate  samples  listed  in  Table  2  are  presented 
in  Fig.  3.  These  values  were  used  only  to  de¬ 
termine  the  pressure  at  which  the  hydrate  sam¬ 
ples  exited  their  stability  regime,  and  thus,  cor¬ 
rection  for  the  elastic  waveguide  effect  of  the 
borosilicate  tube  is  unnecessary.  Dissociation 
pressures  were  determined  by  interpolating  the 
pressure  at  which  cm  fell  below  95%  of  the 
mean  of  the  values  of  cm  corresponding  to  the 
ten  highest-pressure  data  points  for  each  sam¬ 
ple.  Dissociation  was  confirmed  by  a  presence  of 
bubbles  observed  via  the  video  system  at  these 
pressures.  Acoustically-determined  dissociation 
pressures  are  presented  in  Fig.  4  alongside  pre¬ 
dictions  determined  by  several  common  hydrate 
stability  models.  The  Parrish  and  Prausnitz 
model  applies  only  to  structure  I  methane  hy¬ 
drates,  and  thus  does  not  appear  for  comparison 
to  experimental  data  of  samples  CMsII  and  GM- 
sll. 

Data  presented  in  Fig.  3  were  analyzed  for  rela¬ 
tive  changes  in  mixture  sound  speeds,  thus,  cor¬ 
rection  for  the  elasticity  of  the  borosilicate  tube 
was  unnecessary  for  dissociation  pressure  mea¬ 
surements.  Bulk  moduli  measurements  were  de¬ 
termined  using  values  of  cm  obtained  for  sta¬ 


ble  hydrate  samples  (where  no  bubbles  were 
present  in  the  tube),  correcting  cm  to  account 
for  the  elastic  waveguide  effect,  then  calculating 
the  sample  bulk  moduli  using  Eq.  6. 

DISCUSSION 

In  Fig.  4  it  is  shown  that  for  structure  I  hy¬ 
drate  samples,  the  Parrish  and  Prausnitz  pre¬ 
dictions  appear  to  be  outliers  among  the  model 
predictions,  possibly  due  to  the  sample  tempera¬ 
ture,  which  in  the  present  experiment  was  below 
the  0°C  to  27°C  range  validated  in  their  paper. 
Neglecting  the  Parrish  and  Prausnitz  model,  it 
is  evident  that  the  model  predictions  for  sam¬ 
ples  CMsI,  CMsII,  and  HMsI  trend  lower  than 
dissociation  pressures  determined  acoustically. 
This  discrepancy  is  likely  caused  by  the  use  of 
a  strong  brine  in  the  experiment,  as  this  level 
of  salinity  is  not  taken  into  account  in  any  of 
the  models.  The  models  may  have  failed  to  ac¬ 
curately  predict  dissociation  pressures  because 
they  did  not  fully  account  for  the  activity  of 
the  fluid  molecules  surrounding  the  sample.  Al¬ 
though  the  presence  of  simple  salts  does  not  di¬ 
rectly  affect  the  hydrate’s  enthalpy  of  formation, 
dissolved  ions  in  the  fluid  decrease  the  entropy 
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sample 

structure 

/meas  [kHz] 

P  [kg/m3] 

K  [GPa] 

CMsI 

si 

0.1  to  3.5 

1029 

2.8 

CMsII 

sll 

0.1  to  3.5 

2050 

0.4 

HMsI 

si 

0.1  to  3.5 

1497 

0.3 

GMsII 

sll 

0.1  to  3.5 

1326 

1.7 

Lee  et  al. 

si 

500  to  800 

920 

6.4 

Helgerud  et  al. 

si 

not  given 

900 

7.9 

Waite  et  al. 

si 

1000 

900 

7.7  ±0.5 

Pandit  and  King 

si 

820 

750  to  850 

5.7  ±0.1 

Table  3:  Properties  of  methane  hydrate  samples  of  the  present  experiment  compared  to  values  deter¬ 
mined  by  Lee  et  al.  (Ref  [18]),  Helgerud  et  al.  (Ref.  [19]),  Waite  et  al.  (Ref.  [20]),  and  Pandit  and 
King  (Ref.  [21]).  Frequency  range  used  in  measurement  is  indicated  by  /me as- 


of  the  water  molecules,  thus  decreasing  dissocia¬ 
tion  temperatures  [15].  As  temperature  was  held 
constant  in  the  present  experiment,  and  hydrate 
dissociation  pressure  and  temperature  have  an 
inverse  relationship,  the  high-salinity  brine  ef¬ 
fectively  increased  the  pressure  at  which  the  hy¬ 
drate  samples  became  unstable.  A  striking  ob¬ 
servation,  then,  is  each  of  the  models’  seemingly 
accurate  predictions  of  the  dissociation  pressure 
of  sample  GMsII.  The  sample  was  suspended 
in  a  brine  of  132%o  salinity,  yet  all  three  mod¬ 
els  predict  dissociation  pressures  within  6%  of 
the  acoustically-determined  value.  Here,  we  are 
likely  seeing  not  impressive  accuracy  of  the  mod¬ 
els,  but  a  coincidence  of  two  opposing  errors. 
The  first  error  is  caused  by  the  salinity  of  the 
host  solution,  which  raises  dissociation  pressure, 
and  the  second  error  relates  to  the  chemical  com¬ 
position  of  the  gas  in  the  hydrate  samples.  An 
overview  of  the  major  components  of  the  gas  in 
the  samples  is  given  in  Table  1.  While  the  gas 
in  samples  CMsI  and  HMsI  was  almost  purely 
methane,  and  over  80%  of  the  gas  in  sample  CM- 
sll  was  methane,  sample  GMsII  contained  a 
gas  comprised  of  less  than  30%  methane.  Lu 
et  al.  experimentally  found  a  relationship  be¬ 
tween  the  gas  composition  and  the  stability  of  a 
hydrate,  showing  that  stricture  I  hydrates  con¬ 
taining  nearly  pure  methane  were  the  least  sta¬ 
ble,  while  structure  II  hydrates  with  significant 
propane  and  butane  constituents  dissociated  at 
much  lower  pressures  [9]. 

Although  several  models  are  available  to  de¬ 
scribe  the  temperature-pressure  relationship  of 
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hydrate  stability,  the  models  have  limited  ranges 
of  validity,  and  often  do  not  account  for  the 
structure  of  the  hydrate  or  its  surrounding 
medium.  Moreover,  many  of  the  available  mod¬ 
els  have  been  developed  as  fits  to  a  relatively 
small  bank  of  experimental  data.  Further  devel¬ 
opment  of  stability  models  and  new  laboratory 
investigations  of  hydrate  behavior  are  necessary 
to  fully  map  the  phases  of  hydrates  in  arbitrary 
media. 

In  addition  to  measurement  of  hydrate  dissoci¬ 
ation  pressures,  bulk  moduli  of  the  stable  gas 
hydrate  samples  were  measured  using  Wood’s 
model  of  a  two-phase  mixture  (Eq.  6).  Measured 
values  are  presented  alongside  values  found  in 
the  literature  in  Table  3.  Inspection  of  the  data 
shows  values  of  low-frequency  bulk  moduli  which 
are  far  below  previously  measured  values.  Three 
possible  explanations  are  presented  to  account 
for  this  discrepancy.  First,  during  the  experi¬ 
ment  the  video  monitor  used  to  measure  brine 
column  height  became  clouded  with  an  oily  film 
which  was  released  from  the  hydrate  samples. 
This  may  have  resulted  in  inaccurate  measure¬ 
ments  of  hydrate  sample  density  and  volume 
fraction — two  key  terms  in  the  expression  used 
to  calculate  the  bulk  modulus.  For  the  four  sam¬ 
ples  of  the  present  experiment,  an  estimate  of 
the  length  uncertainty  is  ±2  mm  in  liquid  level, 
which  would  result  in  12%  to  28%  error  in  bulk 
modulus  measurement. 

The  oily  substance  released  from  the  hydrate 
samples  is  an  indicator  of  a  second  potential 
source  of  error — impurity.  A  pure  gas  hydrate 
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would  contain  only  H20-ice  and  gas  molecules, 
and  should  have  an  appearance  similar  to  or¬ 
dinary  water  ice,  yet  the  samples  tested  in  the 
present  work  visibly  contained  a  significant  mud 
and  clay  content.  This  could  explain  the  high 
density  measurements  of  the  samples.  However, 
the  literature  shows  that  values  of  sediment  bulk 
moduli  are  typically  similar  to  or  greater  than 
the  established  values  we  have  seen  for  gas  hy¬ 
drates  [22,  23,  24],  Wood’s  model  defines  the 
effective  bulk  modulus  of  a  medium  as  a  lin¬ 
ear  combination  of  the  bulk  moduli  of  the  con¬ 
stituents  of  the  medium.  Therefore,  any  pres¬ 
ence  of  high-bulk-modulus  sediment  in  a  sample 
should  result  in  a  measurement  of  the  effective 
bulk  modulus  of  the  sample  which  is  equal  to  or 
greater  than  the  bulk  modulus  of  the  hydrate  in 
the  sample.  Thus,  we  can  assume  that  the  low 
measured  values  of  bulk  moduli  were  not  a  result 
of  the  presence  of  sediment  in  the  samples. 

The  most  plausible  explanation  for  the  differ¬ 
ences  we  see  between  bulk  moduli  measured  in 
the  present  experiment  and  values  from  the  lit¬ 
erature  stems  from  the  very  motivation  of  this 
work:  It  is  known  that  gassy  sediments  and  com¬ 
posites  can  exhibit  highly  frequency-dependent 
behavior,  such  that  the  acoustic  velocity  of  a 
gassy  sediment  at  high  frequencies  may  be  more 
than  an  order  of  magnitude  greater  than  its 
acoustic  velocity  below  the  resonance  frequency 
of  the  gas  bubbles  in  the  sediment  [25,  26].  From 
the  expression  for  the  speed  of  compressional 
sound  propagation  through  a  solid, 


where  G  is  the  shear  modulus  of  the  solid,  we 
see  that  a  frequency-dependent  sound  speed  im¬ 
plies  highly  frequency-dependent  elastic  mod¬ 
uli.  For  example,  the  tenfold  difference  between 
low-  and  high-frequency  sound  speeds  which 
has  been  observed  in  gassy  sediments  relates 
to  a  hundredfold  difference  in  elastic  moduli. 
Given  the  vast  difference  between  the  frequen¬ 
cies  used  in  the  present  experiment  and  the  mea¬ 
surement  frequencies  found  in  the  literature,  it 
is  plausible  that  all  of  the  bulk  moduli  listed 
in  Table  3 — the  established  high-frequency  val¬ 
ues,  as  well  as  the  low-frequency  values  of  the 


present  experiment — are  correct.  It  is  likely  that 
the  frequency-dependent  behavior  of  the  elas¬ 
tic  moduli  of  methane  hydrates  is  the  primary 
cause  of  the  differences  between  bulk  moduli 
which  have  been  measured  at  ultrasonic  frequen¬ 
cies  and  the  values  measured  in  this  work. 

Additional  Work 

Despite  the  known  frequency  dependence  of  the 
acoustic  properties  of  gassy  substances,  and  de¬ 
spite  a  discrepancy  between  the  low  frequencies 
used  to  locate  gas  hydrates  and  the  high  frequen¬ 
cies  most  often  used  to  characterize  them,  few 
experiments  have  investigated  the  low-frequency 
acoustic  properties  of  gas  hydrates. 

The  promising  results  of  this  experiment  give 
cause  to  continue  work  on  low-frequency  acous¬ 
tics  of  gas  hydrates.  The  authors  intend  to  per¬ 
form  similar  measurements  on  synthetic  si  and 
sll  hydrate  samples.  The  purity  of  synthetic 
samples  will  allow  for  more  accurate  characteri¬ 
zation  of  the  low-frequency  properties  of  the  ma¬ 
terials  and  better  calibration  of  the  measurement 
technique.  In  future  experiments,  temperature 
and  hydrostatic  pressure  will  be  varied  indepen¬ 
dently  to  more  fully  map  the  phase  boundaries  of 
gas  hydrates  using  low-frequency  acoustic  tech¬ 
niques. 
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